Subcritical crack growth: the microscopic origin of Paris's law 
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We investigate the origin of Paris's law, which states that the velocity of a crack at subcritical load 
grows like a power law, da/dt ~ (AK) m , where AK is the stress intensity factor amplitude. Starting 
from a damage accumulation function proportional to (Atr) 7 , Act being the stress amplitude, we 
show analytically that the asymptotic exponent m can be expressed as a piecewise-linear function 
of the exponent 7, namely, m = 6 — 27 for 7 < j c , and m = 7 for 7 > j c , reflecting the existence of 
a critical value 7 C = 2. We performed numerical simulations to confirm this result for finite sizes. 
Finally, we introduce bounded disorder in the breaking thresholds and find that below 7 C disorder 
is relevant, i.e., the exponent m is changed, while above 7 C disorder is irrelevant. 

PACS numbers: 62.20.mm, 46.50.+a, 64.60. av 



In 1963 Paris and Erdogan [lj postulated that under 
fatigue loading a subcritical crack grows with a velocity 
that increases with the stress intensity factor, or equiv- 
alently the crack length, as a power law with an empir- 
ically determined exponent m. Numerous experiments 
showed that this law is valid over at least three orders of 
magnitude for a very wide spectrum of materials The 
Paris law (also known as the Paris-Erdogan law) had huge 
implications in engineering since it allowed to predict the 
residual lifetime of loaded materials quantitatively. To- 
day this law constitutes part of basic knowledge and is 
taught in elementary courses on mechanics jjjj]. Although 
there has been attempts to derive the Paris law in terms 
of geometrical and crack-tip damage-accumulation mod- 
els (see e.g. Ref. [j| and references therein), no work has 
been capable of establishing a firm connection between 
the Paris exponent and microscopic parameters. It is the 
aim of our paper to present analytical and numerical cal- 
culations relating the Paris exponent m to the local dam- 
age accumulation law. This constitutes a micro-macro 
derivation of the Paris law. 

For bodies under cyclic load with stress range Act, sub- 
critical fatigue crack growth follows the Paris law 



da 
~d7 



(AK) m — C (AaVa)' 



(1) 



where a is the crack half-length, AK ~ Aa^fa is 
the stress-intensity factor range, and m is a material- 
dependent exponent. Integration over time leads to the 
Basquin law [J], tf ~ (Act) m , which relates the life- 
time tf (or equivalently the number of loading cycles 
to failure) to the stress amplitude. Let us note that 
while the Basquin law applies to high-cycle fatigue with 
an exponent m depending on the material structure, at 
low-cycle fatigue the corresponding empirical relation, 
tf ~ (Ae)~ A , where Ae is the plastic deformation, is 
called the Coffin-Manson law and has an exponent A 
that is remarkably close to 2, at least for polycrystalline 
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Figure 1: A very thin crack of length 2a, propagating along 
the dashed line in a planar medium subject to a stress 00 at 
infinity. Point P is at a distance x from the crack center. 



single-phased metals Recent work Q has shown that, 
in the rupture of fiber-bundle models subject to fatigue 
damage governed by damage-accumulation functions of 
the form (Ac) 7 , the Basquin law is verified with an ex- 
ponent given by the damage-accumulation exponent 7. 
This brings the question as to whether it is possible to es- 
tablish a direct connection between a microscopic damage 
accumulation form and the Paris law, via the correspond- 
ing exponents. Here, the term 'microscopic' denotes the 
micrometer rather than the atomic (nanometer) length 
scale. 

In order to address this question, we consider a linear 
crack of length 2a, in a two-dimensional medium subject 
to a transverse external stress <To exerted very far from 
the crack, as depicted in Fig. [TJ We model the medium 
as composed of small elements connected by stiff elastic 
springs, with separation Sr. In the continuum limit, and 
within linear elasticity theory, the transverse stress at a 
point along the crack line, a distance x from the midpoint 
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of the crack, is given by 0] 
a {x; a) = <jq 



y/x — a\fx + a 



(2) 



Close to the cr ack tip (x ~ a), the stress diverges as 
a(x;a) = Kj \j2it (x — a), defining the stress intensity 
factor K = aQy/¥a. 

We assume that the medium is under cyclic load, 
with an external stress varying between a m - m and 
Cmax, leading to a stress-intensity-factor range AK = 
(cmax — cr m i n ) yfWa. We further assume that fatigue dam- 
age is the sole factor driving crack growth, which happens 
only along the crack line [8| . Specifically, the half-length 
of the crack increases by Sr when the accumulated dam- 
age at the crack tip reaches a threshold value F t hr- Dam- 
age increments are assumed to be given by 



SF (x; a) = f St (a) [Act (x; a)]" 



(3) 



where St (a) is the number of cycles during which the 
crack has length 2a, /o is a constant related to the time 
scale, and Act (x; a) is calculated from Eq. ([2l with a 
varying between er m i n and £r max - The heuristic [9] power- 
law dependence of the damage increment can be justified 
by invoking concepts of self-similarity 1 1 Oil . and 7 can be 
seen as a stress-amplification exponent [111 ]. 

The accumulated damage at point x when the crack 
has length 2a is given by 



F (x; a) = F(x;a- Sr) + SF (x; a) 



(4) 



Since the half-length of the crack increases from a — Sr 
to a when the accumulated damage F (a + Sr; a) reaches 
Fthr, it follows from Eqs. ([3| and that 



St (a) 



Fhr - F (a + Sr; a - Sr) 
f [Aa(a + 5r;a)F ■ 



(5) 



Here, in order to avoid the appearance of infinities, we 
assume that the spring attached to an element at position 
x experiences a stress given by a (x + Sr;a). To a first 
approximation, this is consistent with the fact that lin- 
ear elasticity theory must break down in the immediate 
vicinity of the crack tip giving rise to a fracture process 
zone or plastic zone [3l.ll2|. 

Numerical iteration of the above equations reveals a 
time dependence of the crack length 2a which, for large 
values of a, reproduces the Paris law, Eq. {TJ), with a 
7-dependent exponent m, as shown in the inset in Fig. 
El Notice that m (as determined from the slopes of the 
log-log plots) seems to reach a minimum value for 7 ~ 
2. This can be explained as a result of the competition 
between damage accumulation mostly around the crack 
tip, which happens for 7> 1, and uniform damage along 
the whole crack line, which dominates as 7 approaches 
zero. The existence of this minimum is confirmed by 
calculations of m as a function of 7, obtained by power- 
law fits of da/dt versus a, using various system sizes L. As 




Figure 2: (Color online) Inset: Log-log plot of the time deriva- 
tive of the crack half-length o, as a function of a/ao, for var- 
ious values of the damage exponent 7. In each case, the time 
scale was adjusted so that da/dt approaches unity as a ap- 
proaches the initial value ao = 100<5r. Main panel: Depen- 
dence of the Paris exponent m on 7, for different values of the 
system size L. The dashed curve corresponds to Eq. ([6]). 
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Figure 3: (Color online) Scaling plot of m versus 7 for different 
values of the system size L, showing conformance with Eq. 
([7]). The best data collapse is obtained with x = 0.089. 



shown in the main panel in Fig. [21 there are strong finite- 
size effects around j c = 2. As L — > 00, the minimum 
in each curve shifts slowly towards m = 2 at 7 C = 2 , 
suggesting an asymptotic form 



6 — 27, for 7 < 7c 
7, for 7 > 7c 



(6) 



This can be checked by a finite-size scaling hypothesis, 
1 



m (L) 



L*F+(|7- 



2\L X ) 
2\L X ) 



for 7 < 7c 
for 7 > 7c 



(7) 



which, as shown in Fig. [31 is nicely fulfilled by our numer- 
ical data with x — 0.089. Both scaling functions behave 
as T± (u) r~j u^ 1 , for u ^> 1, in agreement with the sug- 
gestion that, in the continuum limit [L ^> Sr), m should 
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be a linear function of 7 for both 7 < j c and 7 > j c . We 
have therefore found evidence that 7 = 7 C is a critical 
point, associated with the divergence of the stress inte- 
gral along the crack line. Our conclusion that m > 2 is 
fully compatible with values reported for various materi- 
als throughout the literature. 

The prediction of Eq. is confirmed by an analyt- 
ical treatment of Eqs. ([H])-©. From the time St(ao) 
during which the crack length is 2ao, we can calculate 
the instantaneous damage F (ao + 2Sr; 00) at the next 
crack tip position eto + 2Sr, and then recursively deter- 
mine the damage fraction at the successive crack tip posi- 
tions, G n = F (a + Sr; a — Sr) /F t hi, from St (a — Sr) and 
F (a + Sr; a — 25r), with a — a Q + nSr. As a result, it is 
possible to write the recurrence relation 



fc=i 



fc=i 



9n,k, 



with Go = and 

'Act (oo + (n + l)Sr; a + (fc - l)Sr) 



9n,k 



1 T 



Act (ao + kSr; ao + (k — l)Sr) 



(8) 



(9) 



Notice that is determined by the ratio between the 
stress amplitudes at position nSr and at the crack tip, 
when the crack has grown by a distance 2kSr. From Eqs. 
([5]) and ([2]), it follows that, for a = ciQ + n5r, with 



da 



Sr 
St (a) 



1-G n 



(10) 



so that the scaling behavior of the crack growth rate de- 
pends on the scaling behavior of G„. This can be inves- 
tigated by looking at Eqs. J8|) and ([9]), from which, if 
Sr <C ao -C nSr, we obtain the asymptotic forms 



9n,k 



(2Sr/a )^ : 
(n-k + 2Y 



-1 



for k <C ao/Sr; 

for ao/Sr <C k <C n; 

for fc < n. 



(11) 



It is also easy to show that g n ,k has a single minimum 
at fcmin ~ n/y/3, so that, as n — > 00, the sum on the 
right-hand side of Eq. ([8]) has a power-law divergence 
n 1_ 2T for 7 < 2, while it converges to a finite value 
C 7 = C(It) — 1 f° r 7 > 2, where £(x) is the Riemann zeta 
function. It follows that G n approaches C 7 /(l + C 7 ) < 1 
for 7 > 2, and from Eq. iflOj) we immediately see that 
rn = 7. On the other hand, numerical calculations show 
that, for 7 < 2, G n asymptotically approaches unity; al- 
though we were not able to derive an analytic expression, 
it can be readily checked numerically that in this case 



1-G„ 



-3(i-H 



(12) 



leading, when combined with Eq. (fT0|) . to m = 6 — 27. 

For 7 > 7 C = 2, in the light of Basquin's law, the 
prediction m — 7 is compatible with the result obtained 



in Ref. [6j], stating that the Basquin-law exponent is 
given by the damage exponent 7. This is no longer valid 
for 7 < 7 C . However, the models studied in Ref. 
involve randomness in both fatigue and stress thresholds 
as additional ingredients. 

In order to investigate the effects of disorder on the 
Paris exponent, we introduce a distribution of values of 
the fatigue thresholds with lower (upper) cutoff Fa (F2), 
so that each lattice point has a local threshold Fi < 
Fthr(x) < F2. An obvious consequence of the disorder 
is that, depending on the disorder strength and on the 
damage exponent 7, it is possible that a point far from 
the crack tip reaches its local fatigue threshold before a 
point closer to the crack tip. This leads to the occurrence 
of rupture avalanches, in much the same way as in fiber- 
bundle models (see e.g. Ref. [13]). 

Conditions for the appearance of avalanches can be 
derived from the analytical approach presented above. 
Avalanches will occur if, for some crack length 2a, the 
local threshold at position a + Sr is less than the accumu- 
lated damage at that position when the crack had length 
2(a - Sr), i.e. F thr (a + Sr) < F (a + Sr; a - Sr). As the 
fatigue thresholds are no longer the same for all points, 
Eq. J5]) ceases to be valid. However, F (a + Sr; a — Sr) 
can still be written as a linear combination of the thresh- 
olds of the points closer to the crack center, with coeffi- 
cients related to the stress ratios g n ,k- In the absence of 
any previous avalanches, and for 7 > 7 C = 2, the value of 
F (a + Sr; a — Sr) is limited by F2 times the asymptotic 
value of G„, i.e. F (a + Sr; a - Sr) < F 2 C 1 j (1 + C 7 ). 
Since F\ < F t hr (a + Sr) , we conclude that avalanches 
will occur if 



F) 



-T2 



1 + C7 7 



< 



F\ < C 7 
F 2 ~ 1 



Cry 



(13) 



On the other hand, for 7 < 7 C , any finite amount of 
disorder leads to the occurrence of avalanches. 

Confirmation of this prediction, as well as further in- 
formation on the effects of disorder, can be obtained from 
numerical calculations. In order to analyze the results of 
those calculations, we integrate the Paris law to obtain 



1 - (a /a) 1 



Bt, 



(14) 



where B is a m-dependent constant related to the inverse 
rupture time. Introducing fatigue thresholds uniformly 
distributed between F\ and F%, numerical iteration of 
Eqs. ((2|) and ^ — modified to acomodate local thresh- 
olds — shows that, for 7 > j c , plots of 1 — (ao/a)^ m_1 
for large times (not shown) remain straight lines, with 
the same value of m = 7 as in the absence of disorder; 
however, for sufficiently small values of / = F1/F2, the 
coefficient B becomes disorder-dependent. 

For 7 < 7c, we confirm the occurrence of avalanches for 
any amount of disorder. Also, as shown in Fig. |4] (main 
panel), the prediction of Eq. lfT4|) is no longer verified, 
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Figure 4: (Color online) Rescaled plots of the crack length for 
7 = | and various ratios / = F1/F2 between the minimum 
and maximum values of the fatigue thresholds. Main panel: 
rescaling as given by Eq. lfl4|l . Inset: rescaling with an expo- 
nent 1.14 7^ 2 m ~ 1- AH curves were obtained for system size 
L = 2 15 and averaged over 500 disorder realizations. 

indicating that the Paris exponent m is actually modified 
by the introduction of disorder The new exponent 
to', which recovers the linear behavior predicted by Eq. 
(fl4|) . is almost independent of the ratio /, for sufficiently 
strong disorder, as also shown in Fig. [4] (inset), and we 
have checked that in this limit the ratio {ml — 2)/(m — 2) 
is independent of 7, being given approximately by 0.76. 

We have been able to derive analytically and confirm 
numerically the Paris law and found that its exponent m 
is a function of the damage exponent 7 describing the mi- 
croscopic damage accumulation law. To our surprise we 
discovered that j c = 2 is a critical point characterized by 
a scaling law and a critical exponent which separates two 
regimes with different linear functions to (7). In addition 
we also studied the role of disorder and found again that 
7 C = 2 plays a special role: disorder is relevant below it 
and irrelevant above it. Our results can have far-reaching 
consequences in the understanding and control of subcrit- 
ical crack propagation. On one hand the discovered rela- 
tion between the damage and the Paris exponents, which 
in principle could be checked experimentally, could help 
to predict lifetimes of samples by studying the velocity of 
small cracks. As we found, the value j c = 2 is very spe- 
cial. For example, to avoid the influence of disorder one 
must try to stay above j c . Finally, it is important to men- 
tion that the exponent 7 of the damage accumulation law 
is material dependent and could therefore play a central 
role in the engineering design to increase the robustness 
and optimize the mechanical performance of the system. 
The derivation of power-law damage-accumulation func- 
tions from atomistic processes would be of much interest, 
and represents a big challenge for future work. 
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